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Abstract 

We propose that the Universe is filled with a massive vector field, non-minimally coupled to 
gravitation. The field equations of the model are consistently derived and their application to 
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■ cosmology is considered. The Friedmann equations acquire an extra dark-energy component, which 
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is proportional to the mass of the vector particle. This leads to a late-time accelerated de Sitter 
type expansion. The free parameters of the model (gravitational coupling constants and initial 
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value of the cosmological vector field) can be estimated by using the PPN solar system constraints. 
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I. INTRODUCTION 



An increasing number of observational data, obtained in the past few years, strongly 
support a model according to which the Universe is spatially flat, mostly made of non- 
conventional matter - baryons being allowed only up to 4% of the total energy content 
- and accelerating. The physical models accounting for such a picture generally contain 
two basic ingredients: pressureless dark matter (DM), responsible for the growth of cosmo- 
logical perturbations via gravitational instability, and negative pressure dark energy (DE), 
responsible for the accelerated expansion (for recent reviews on the dark energy problem 
see 

The simplest model along these lines is ACDM, in which the role of DE is played by 
a cosmological constant A. It fits very well all the data related with the cosmological 
background and the perturbations in the linear regime 2j. The cosmological constant is 
attributed to the quantum zero-point energy of the particle physics vacuum, with a constant 
energy density p, pressure p and an equation of state w = p/ p = — 1. Existing observational 
data indicate the equation of state of dark energy very close to the cosmological constant 
value, w = —1 ± 0.2 at 95% confidence level, with at most a very mild evolution up to 
redshift z ~ 1. 

An alternative model to the cosmological constant may be the quintessence a dynam- 
ical scalar field which, in the simplest model, slowly rolls in a potential characterized by 
an extremely low mass. The spatially averaged equation of state for the quintessence field 
satisfies w > — 1. In order to avoid fine-tuning on the initial conditions, the quintessence 
scalar field is usually taken to be extremely light, with a Compton wavelength corresponding 
to the present value of the Hubble radius. As a consequence, the scalar field is homogeneous 
on all observable scales, much like a cosmological constant. However, quintessence models 
suffer from an important theoretical problem, namely, the fact that radiative corrections 
induced by the couplings with the matter fields would generically induce huge corrections 
to the trace level mass, thus spoiling the required lightness. Hence to keep the scalar field 
light in these models a fine-tuning on the radiative corrections is generally required, besides 
the one necessary to keep the cosmological constant small. 

A wide variety of other dark energy models has also been proposed, including K- 
essence (4|, Chaplygin gas [5], modifications of gravity [6], Born-Infeld scalars (rolling 
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tachyon) massive scalars js| etc. The common feature of these models is that they 
operate through an undetermined field potential which in principle can incorporate any a 
priori associated cosmological evolution, thus lacking predictive power at the fundamental 
level tJ. There is a tremendous degeneracy in these models and generally they are judged 
by their physical implications and by the generic features which arise in them. Therefore 
a consistent physical picture of the dark energy, which could explain the size of its energy 
density pde ~ 10~ 12 eV 4 , and suggests how the underlying physics may be probed, is still 
missing. 

Non-gravitational interactions are known to be mediated by vector fields. Therefore the 
possibility that a vector field, which, for example, may be a partner of quintessence, could be 
at the origin of the present stage of the cosmic acceleration cannot be neglected. Theoretical 
proposals in which a minimally coupled vector field is responsible for the present dynamics 
of the Universe have been considered in 9|. 

It is the purpose of the present paper to propose a model of the dark energy in terms 
of a massive, Proca type vector field, with a non-minimal coupling to the gravitational 
field. The model contains three independent parameters u, r\ and fi\, respectively, with 
n\ representing the mass of the massive cosmological vector particle. The scalars to and r\ 
describe the non-minimal coupling of the vector field to the Ricci scalar and to the Ricci 
tensor, respectively. The gravitational field equations can be consistently derived from 
a variational principle. A similar vector-tensor theory, without the mass term, was also 
proposed in the early 1970's flo| . 

In the cosmological case, corresponding to a flat homogeneous and isotropic Universe, 
the Friedmann equations acquire an extra dark-energy component, which is proportional to 
the mass of the vector particle. This term, playing the role of the cosmological constant, 
leads to the late accelerated expansion of the Universe. The free parameters of the model 
(gravitational coupling constants and initial value of the cosmological vector field) can be 
estimated by using the PPN solar system constraints. The mass of the cosmological massive 
vector particle, which may represent the main component of the Universe, is of the order of 

1Q -63 g 

The present paper is organized as follows. The field equations of the massive vector-tensor 
theory are derived in Section II. In Section III we consider the cosmological applications of 
the model. The PPN constraints on the model parameters are discussed in Section IV. We 
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discuss and conclude our results in Section V. 

Throughout this paper we use the Landau-Lifshitz conventions 11( for the metric signa- 
ture (+, — , — , — ) and for the field equations, and a system of units with c = h = 1. 



II. FIELD EQUATIONS OF THE MASSIVE VECTOR-TENSOR THEORY 

We assume that the Universe is filled with a massive cosmological vector field, with mass 
/ia, which is characterized by a four-potential A IJ, {x u ), v — 0,1,2,3 and which couples 
non-minimally to gravity. In analogy with electrodynamics we introduce the field tensor 



(1) 



The interaction of the gravitational and of the vector fields is described by a Lagrangian 
which is required to satisfy the following conditions: a) the Lagrangian density is a four- 
scalar b) the free-field energies are positive-definite for both the metric and the vector field 
c) the resulting theory is metric and d) the field equations contain no higher than second 
derivatives of the fields jlOj] . The action for such a theory can be written as 

1 

2' 



S 



R + C^W + ^/i A A^ + cuA^R + rjA^R^ + 16nG L r 



-gdn, (2) 



where R^ v and R are the Ricci tensor and the Ricci scalar, respectively, Go is the gravi- 
tational constant and L m is the matter Lagrangian. In Eq. (j2J) uj and r\ are dimensionless 
coupling parameters. The four- dimensional volume element is dQ = dx°dx 1 dx 2 dx 3 . In the 
following we denote <fi = A^A^, which is an invariant scalar. 

The variation of the action with respect to the metric tensor g^ v gives the field equations 



R 



-g^R + uj 



R 



-g^R 1 + A„A,.i?- (/ /( ,V X V ; 



V„V, 



r]A a A^ ( 9lm R 



[g,uV a V p {A a A?) + V a V° (A M A„) - V^V, (A M A°") - V a V, (A„A<0] + 
ZCfiaCu ~ -jg^CupC^ + -/i^A^A;, — ~n\(f>g[i,v = SttGqT^, 



(3) 



where T^ v is the energy-momentum tensor of the matter, defined in terms of the matter 
action S m = 167rG / L my /—gdD, as SS m = 8ttG J T^Sg^^—gdft. 
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The variation of the action with respect to A M gives the generalized Maxwell equation 

AS7 U C» V + uRA^ 1 + rjA"R£ + -^A" = 0. (4) 

By contracting the field equations we find 

- R + (3u + |) V A V> + r/V a V/j (A a A^) = 8ttGT, (5) 

where T = T^. Taking the covariant derivative of Eq. (jlj) we obtain the conservation law 
for the four-potential of the massive cosmological field as 

Q/4 + uR) V m A" + (lu + |) (V M i2) A" + r/ (V M A") Jtf = 0. (6) 

From Eq. (jlj) it follows that the four-potential of the cosmological vector field satisfies 
the following wave equation: 



V a V a A» - V M (VA V ) - 



A" = 0. (7) 



Due to its antisymmetry, the massive vector field tensor automatically satisfies the equa- 
tions 

VaC^ + V^C va + V v C oll = 0. (8) 

The matter energy momentum tensor satisfies the conservation law V^T^ = 0, which 
can be verified by taking the covariant divergence of Eq. 

III. COSMOLOGICAL APPLICATIONS 

To investigate the cosmological implications of the massive vector-metric theory we adopt 
the flat Robertson- Walker metric for a homogeneous and isotropic Universe, given by 

ds 2 = dt 2 - a 2 it) (dx 2 + dy 2 + dz 2 ) , (9) 

where a(t) is the scale factor. The observed isotropy and homogeneity of the Universe 
requires that the massive vector field is a function of the cosmological time only. Hence we 
assume that the potential A M has only one non-zero component, A M = (A°(£), 0, 0, 0). The 
function <fi is given by </>(t) = A (t)A°(£). 

The non-zero components of the Ricci tensor and the Ricci scalar are given by Rqo = 
— 3d/a, R aa = (ad + 2d 2 ) S aa , a = 1, 2, 3 and R = — 6 (a/a + a 2 /a 2 ), respectively. Moreover, 



we have V a V/3 (A a A^) = + 30a/a + 60d 2 /a 2 + 60d/a, V^V 7 (A A ) = + 30d/a-60d 2 /a 2 
and V^Vo (A A CT ) = 4> + 3<pa/a - 30d 2 /« 2 - 

We assume that the matter content of the Universe consists of the massive cosmological 
vector field C^ u and ordinary matter in form of pressureless dust with density p. The 
conservation of the energy- momentum tensor then gives p = po/a 3 , where po is a constant 
of integration. 

Hence, the gravitational field equations and the equation of motion of the vector field 
become 

r-i / \ n a 2 v. , Oj ( 77 \ -a u\. 8nG p ... 

[2 + (2c + 3^)0] ^ + [l + (^ + 37^)0]^+^ + + 3(^ + 7/)^ -^0 = 0, (11) 

a a 2 6 

Eliminating with the use of Eq. (fl2l) the second derivative of the scale factor from Eq. (flUj) 
gives 



[l + (u, + ,)0]-+ (u> + -j = — - (13) 

The general solution of Eq. (fl2l) can be represented in an integral form as 



where to and ao are arbitrary constants of integration. 

In the limit of large a and for Auj (2u + if) > the scale factor is given by 

a = exp [H (t — t )] , H = constant, (15) 

corresponding to a de Sitter type exponentially accelerating phase for the expansion of the 
Universe. The constant H is expressed in terms of the mass of the cosmological vector field 
as 

H 2 _ Pa _ A (IQ) 
6{Auj + t]) 3' { ' 

where A is the cosmological constant, which is generated due to the presence of the massive 
cosmological vector particle, and whose numerical value can be determined from observa- 
tions. 

For arbitrary times the exact form of the scale factor is given by 

a(t) =a smh n [(3(t-t )}, (17) 



where we have denoted 



and 



n 



2u + 7] 



Au + 7] 

(4w + 7])H H 



(18) 
(19) 



(2oj + 7}) n ' 
respectively. 

With the use of Eq. (|T7|) it follows that the massive cosmological vector field satisfies the 
evolution equation 

2H 



4 = p 



tanh (3(t- t ) - 2 p + ^ coth /? (t - t ) 

AuJ + 7] 



coth /3 (t - t ) 



+ 



2u; + 7/ 

167rG po sink 1 - 3 " /3 (t - t ) 
3 (2cu + 77) apf cosh /? (t - t ) 



(20) 



with the general solution given by 



cosh 



\P (t ~ *o)] f 



5 + 



167rG po 



sinh^ [/3 (t _ to)] I 3 (4a; + „) a§IP 



tanh [0 (t - t )] 



AiU + 7] 



F [/?(*- to)]}, (21) 



where .B is an arbitrary constant of integration, and 

dx 



F(x) 



where 



m 



sinh m (x) 

2u — 7] 



(22) 



(23) 



4uj + 77 

During the pure de Sitter phase, with the effect of the ordinary matter neglected, the 
equation describing the dynamics of the time varying cosmological vector field is 



2u — 77 

2UJ + T] 



H<f> 



2u + 7] 



■H, 



with the general solution 

(t) = 0o exp 



2uj — 77 



Ht + 



2uj — 7] 



1 — exp 



2uj — 7] 
2uj + 7] 



Ht 



2u + 77 

where 0o is the initial value of the field at the initial time t — to — 0, (0) - 
In the case of a constant field = 0o =constant, H and are related by 



(24) 



(25) 



H 1 



Ma 



12 1 + (u; + 77) 0o 

from which we find the constant massive cosmological vector field as 

4A//4 2 

^° 1-4(0; +77) A/74 2CJ-77' 



(26) 



(27) 
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IV. PPN CONSTRAINTS ON THE MODEL PARAMETERS 



The numerical values of the coupling coefficients 00 and 77 and the initial value </>o of 
the cosmological vector field can be constrained by using solar system observations. Vector 
tensor models generate observable effects in light deflection and retardation experiments, 
planetary perihelion advance, orbiting gyroscope precession, non-secular terms in planetary 
and satellite orbits, geophysical phenomena etc. These effects can be described in terms 
of the dimensionless parameters a, (3 and 7, which parameterize deviations with respect to 
standard general relativity. Actually, the parameter a can be settled to the unity due to the 
mass definition of the system itself [121 ] . 

The quantity 7 — 1 measures the degree to which gravity is not a purely geometric 
effect, and it is affected by other fields. Measurements of the frequency shift of the radio 
photons to and from the Cassini spacecraft as they passed near the Sun give the result 
7 = 1 + (2.1 ± 2.3) x 10~ 5 [13]. The value of the parameter 4/3 — 7 — 3 can be constrained 
from Lunar Laser Ranging, with the observational result 4/5 — 7 — 3 = — (0.7 ± 1) x 1CT 3 14]. 
For the cosmological time variation of the effective gravitational constant we adopt the 



value G/G = 10 14 yr 1 15]. As for the present value of the Hubble constant we take 
H = 70km/s/Mpc. 

Therefore, the three free parameters of our model (<f) ,uj,ri) can be obtained from the 
following non-linear system of algebraic equations: 

7 -l = ^P^ = 2,lxl0- (28) 



4(3 - 7 - 3 



(7 ~ 1) (1 

U)(j) 



LUC 



1 + 



7 (7 - 2) 



1(7 + 1) 



5(7-1) 



■luJ 



+ 



3^(7-1) I 5(7-3) 



+ 



-0.7 x 10" 



(29) 



1 G 
HG 



3^(7-1) 1 ^(7-3) 
2 ~ r 4 



[(2u;-77)0o-2]exp'^ 



2u + V {!(7 + l)-^i) + 



3^(7-1) 1 ^(7-3) 
2 ' 4 



1.40 x 10" 



where we denoted by </> the present day value of the cosmological vector field, 



(30) 



2u — i] 



exp 



2u — 1] 
2oj + i] J ' 2oj — 77 



(31) 
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These three equations are highly non-linear, and approximate solutions can only be ob- 
tained by the means of numerical methods. 

To constrain the numerical values for the three free parameter, we firstly consider the 
(u, o )-plane for given values of the parameter rj in Figures [JJ and [2J 



?o (w,0 o ) plane for /7 = 18000 



1.5 



0.5 



8990 9000 9010 



'o (OJ,0 O ) plane for /7 = 16000 



1.5 



0.5 



FIG. 1: These figures represent the (u, (/>o)-planes f° r V = 18000 and rj = 16000. The three straight 
lines in the middle represent the allowed range of the parameter 7 — 1 = (2.1 ± 2.3) x 10 -5 , as given 
by Eq. (|28p . The dashed lines represent the allowed range of 4/3 — 7 — 3 via Eq. (|29p . The right 
dashed lines correspond to the value 4/3 — 7 — 3 = —(0.7 — 1) X 10~ 3 = +0.3 x 10~ 3 and implies 
that the model strongly favors positive values. Consequently, a further decrease on the error bars 
of 4/3 — 7 — 3 can either rule out or further strengthen the present model. Lastly, the dotted 
lines describe the allowed parameter range for the 'cosmological' equation (|30p . where we used 
H = 70 ± lOkm/s/Mpc. It is evident that there is no parameter range where all three equations 
are satisfied to given accuracy, in fact one has to search for parameters that minimize the error. 



The (rj, 0o)-pla n es share features very similar to those of the (u, 0o)-planes, therefore 
it is not particularly insightful to also consider these. The reason behind this is that the 
numerical fitting of the three parameters works best if 2u; — 77 is of the order O(10 -5 ). The 
smallness of this difference has its roots in the Newtonian limit of the theory, which in case of 
a massless vector field reduces to 2u = rj. Hence, this small deviation is due to the massive 
vector field. However, let us consider the (u, ??)-planes for given <ft in Figure El 

For example, one possible set of numerical values for the free parameters of the model is 
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(w,0o) plane for (7 = 8000 *o (w,0o) plane for r/ = 1000 



1 1 ]L_L 



3980 3990 4000 4010 4020 



"T rl'i 



495 500 505 510 



FIG. 2: Using the same coding as in the previous figure, it should be noted that the possibility of 
fitting the three parameters becomes more difficult. Already for r/ = 8000 the 'cosmological' range 
of parameters (dotted) lies outside the allowed range of parameter as implied by Eq. (f29j) (dashed). 
Hence, it is even more difficult to fit the three parameters simultaneously. A further increase of r/ 
larger than 18000 has the same effect, the situation does not improve. 



H (w,/7) plane for O =1.OO *7 plane for O =1.OO 




2000 4000 6000 8000 10000 8800 8850 8900 8950 9000 9050 9100 

FIG. 3: As in the above figures, straight lines show the 7 — 1 range, dashed lines the 4/3 — 7 — 3 
range, and dotted lines the cosmological constraint. On the left-hand side, where < uj < 10000, 
the different lines cannot be distinguished. The right-hand side (8750 < lo < 9100) again indicates 
that one can find some approximate solution for the three parameters. An increase of the parameter 
0o does not change the qualitative picture. 
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given by O = 1.097999982, u = 9000.000069565 and T] = 18000.0001381, respectively. These 
values satisfy the general constraints on the coupling constants 2u — r\ > and 4u + r\ > 0. 

The first condition is implied by the positivity of <po, see Eq. ( 12 7p . since is the square 
of the massive vector field. The second condition follows from Eq. ffTBT) . i.e. positivity of 
the vector field mass and positivity of the Hubble constant. The above values represent an 
approximate solution of the system of constraint equations of the order of O (10~ 5 ), O (10 -4 ) 
and O (10~ 4 ). It is interesting to note that the system of equation is quite sensitive towards 
a change of the 'cosmological' equation (|30|) . 

If, for the moment, we assume that the gravitational constant G changes slower than given 
by the upper bound G/G = 10~ 14 yr~ 4 , for example an order of magnitude slower, then also 
the fit of the 'cosmological' equation improves by an order of magnitude. Hence, we can 
conclude that the analyzed theoretical model makes mild predictions regarding observation: 

The gravitational constant should vary slower over time and the parameter 4/3 — 7 — 3, 
constrained by Lunar Laser Ranging, should lie in the positive region. On the other hand, 
it should be emphasized that our model would not be in good agreement with G = 0, a 
constant effective gravitational "constant". 

V. DISCUSSIONS AND FINAL REMARKS 

From Eq. ( fT6l) we obtain for the mass of the cosmological vector particle the expression 



/i A = V6 (4cj + r])H « 1.67 x ^6 {Aw + 77) x 1(T W g. (32) 

The upper limit for the mass of the ordinary photon obtained by using a rotating torsion 
balance method is 1.2 x 10~ 51 g If]]. It should be noted that the proposed massive vector 
field interacts only with gravity and has no standard matter interactions. The existence of 
a minimal mass in nature in the presence of a cosmological constant has been discussed for 

n 

example in [17| . 

The mass of the cosmological massive vector particle may have been generated during 
inflation. Inflation allows the emergence of fields coherent over large distances from quantum 
fluctuations and prevents dissipative effects, due to the absence of charged particles. The 
conformal invariance of the U(l) gauge theory for electromagnetism prevents the gravita- 
tional field from producing photons. However, in the present model, due to the breaking 
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of the conformal invariance, via the coupling to gravity, the vector particle interacting with 
gravity can acquire an effective mass. Therefore, the effective mass for the massive vec- 
tor particle composing the cosmological gas filling the Universe, and representing its main 
matter component, may have originated during inflation }18| . The mass of the cosmological 
vector particle can also be acquired from the spontaneous breaking of Lorentz symmetry in 
the context of field theories arising from string field theory, or due to the introduction of a 



19] . 



fundamental minimal length in a trans-Planckian physics scenario 

The vector field dominated Universe enters in a pure exponential de Sitter phase when 
the condition t > t acc = 1/(3 — (2a; + rj)tH / (4u + rf) is satisfied. With the use of the 
values of the parameters u and rj obtained from the PPN solar system constraints we obtain 
t acc ~ 2t H /3, corresponding to a redshift of z acc « 0. In fact, the supernova data published 
by the High-z Supernova Search Team and the Supernova Cosmology Project show that 
the transition from the decelerating to the accelerating phase occurred at redshifts smaller 
than z = 0.4 j^. Therefore, the model predicts that the accelerated expansion phase of the 
Universe started only recently. 

In conclusion, we have shown that the dark energy can be modeled as a massive cosmo- 
logical vector field filling the Universe. This field may have originated during the inflationary 
period, when the vector field may have acquired its mass, and it drives the late-time acceler- 
ation of the Universe. All the parameters of the model can be constrained from observational 
data. 
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